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The recent interest in the use of pulse code modulation (PCM) for the 
transmission of speech has made it desirable to determine the effects of 
digital line errors on certain codes. Improved performance for low-level 
talkers has been generally obtained by the use of instantaneous compressors 
and expandors (compandors) in order to avoid using a large number of 
digits in the PCM coder-decoder (codec). A comparison of the effects of 
digital errors for a logarithmically compandored system transmitting 
speech is made on the basis of mean square distortion power and the distri- 
bution of error magnitude. Results are obtained for the binary, Gray (re- 
flected binary), and folded binary codes. The results indicate that the binary 
and Gray codes have a "click" as well as a ll noise" component of distortion 
while the folded binary code produces only noise at low-talker levels. "Clicks" 
have been defined as errors which have amplitudes greater than half the full 
range amplitude; the remaining errors are considered "noise". For the 
folded binary code, the most significant digit gives polarity information; 
the remaining digits represent the signal magnitude in binary code. 

I. INTRODUCTION 

A comparison of the effects of digital errors for a logarithmically com- 
pandored 12 system transmitting speech by PCM is made on the basis 
of mean square error power and the distribution of the error magnitude. 
A block diagram of the system considered in the analysis is shown in 
Fig. 1. The calculations are based on the assumption of at most one 
error per 8-digit word. The results obtained for binary and Gray codes 
indicate that there is a "click" as well as "noise" component of distor- 
tion at low talker levels. "Clicks" have been denned as errors which 
have amplitudes greater than half the full range, the remaining error 
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Fig. 1 — A block diagram of the system under consideration. 

amplitudes are considered "noise". A folded binary code (see Fig. 2) 
essentially produces only the noise component at low talker levels. 
Low talker levels were considered to be 50 db below full load sine 
wave (50 db BFLSW) for all the cases calculated. 

The results obtained for the folded binary code indicate that this 
code does not have a click problem at low talker levels when compared 
to the binary and the Gray codes. A question is raised, however, concern- 
ing the subjective effect at low levels of the noise components of dis- 
tortion. If the noise is subjectively equivalent to quantizing error, an 
error rate of 10 -6 will cause a 3-db impairment in the system perform- 
ance. This occurs for speaker levels between 30 and 50 db BFLSW in a 
system using a logarithmic compression factor 1,2 (n) of 100 and the 
folded binary code. 

The conclusions obtained from the computed results for logarithmic 
companding are: 

(i) The folded binary code produces fewer clicks than either the Gray 
or binaiy code at low speaker levels. 

(ii) The Gray, binaiy, and folded binary click and noise components 
are comparable at higher speaker levels. 

The combined subjective effect of noise and clicks occurring simul- 
taneously has not been evaluated. 



II. RESULTS 

The results obtained for the probability of a given code word occurring 
are plotted in Fig. 3 for logarithmic compression factors 1,2 of 50, 100, 
and 200. The input signal was assumed to be Laplace distributed. 3,4,5 
The figure indicates that compression has a major effect on the prob- 
ability of a given code word occurring. The tendency of the compressor 
is to assure a more uniform distribution of the signal across the lower 
code levels for low and average level talkers. One notes that for a /x of 
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Fig. 2 — 5-digit codes. 



100 the average talker (30 db BFLSW) has a nearly uniform distribution 
over the range — £ to +f full range. 

The mean square value of the distortion caused by digital line errors 
is compared to 8-digit quantizing error power 12 for a 10 -6 error rate, a 
2-volt peak-to-peak input amplitude and compression factors (m) of 
50, 100, and 200 in Figs. 4 and 5. Altering the error rate simply results 
in a scale change for the line error distortion power. The results plotted 
in Fig. 4 include all error amplitudes whereas errors greater than half 
the full range amplitude (clicks) have been removed for the n = 100 
case in Fig. 5. Preliminary results of subjective tests indicate that clicks 
are less objectionable than quantizing noise. The figures indicate the 
folded binary code yields superior click performance at low talker levels 
compared to the Gray or binary codes. The noise or non-click component 
for the folded binary is greater than that produced by the Gray code at 
low levels. For an error rate of 10" 6 and 8-digit quantizing, the line error 
distortion power and quantizing error power are equal for the folded 
binary code in the 30 db to 50 db BFLSW range. If the line error dis- 
tortion is subjectively equivalent to quantizing error, the system per- 
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Fig. 4 — Mean square error vs speaker level-8 digit coder — 2-volt peak-to-peak 
full load amplitude, (a) fi = 50 log companding, (b) m = 100 log companding, (c) 
H — 200 log companding. 

formance will be degraded by 3 db. This result applies for a logarithmic 
compression of 100. 

We conclude that the essential low-level problem is noise for the 
folded binary and clicks for the Gray and binary codes. 

The probability of errors of a specified magnitude occurring are 
plotted in Figs. 6 through 9 for the binary, folded binary, and Gray codes. 
A clear demarcation between errors greater than and less than half the 
full amplitude is indicated in Figs. 6, 7, and 8 for low levels. These re- 
sults led to the definition that errors greater than half the full load would 
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Fig. 5 — 8-digit codes (compression factor, /x = 100 error rate 10"°, clicks re- 
moved — 2-volt peak-to-peak full load amplitude). 



be called clicks. The dividing line between clicks and noise is not clear 
for average and high-level speakers. 

It is apparent that the folded binary code produces virtually no 
clicks at low levels. This leads to the conclusion that noise rather than 
clicks will be the dominant low-level problem in transmission of folded 
binary, whereas clicks are the problem for Gray and binary. 

The mean square error caused by specified digit errors for various 
speaker levels are plotted in Fig. 10(a), (b) and (c). The results indicate 
that the largest mean square error for the folded binary code is caused 
by digit 2 errors. The largest mean square error for Gray and straight 
binary depends to some extent on the speaker level. The largest error 
for low and average levels occurs in the 1st digit for binary and the 2nd 
digit for Gray. 
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Fig. 6 — Probability of error (speaker level 50 db BFLSW, fi = 50 log 
companding, symmetric for negative amplitudes, low-level talker), (a) 8-digit 
folded binary code, (b) 8-digit binary code, (c) 8-digit Gray code. 
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Fig. 7 — Probability of error (speaker level 50 db BFLSW, m = 100 log 
compandor, symmetric for negative amplitudes, low-level talker), (a) 8-digit 
folded binary code, (b) 8-digit binary code, (c) 8-digit Gray code. 

III. ASSUMPTIONS AND PROCEDURES 



3.1 Assumptions 

The results obtained were calculated using the following assumptions: 
(i) The terminals of the transmission system including the coder and 
decoder are ideal. 

(it) All digital errors are introduced in the transmission medium. 
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Fig. 8 — Probability of error (speaker level 50 db BFLSW, m = 200 log 
companding, symmetric for negative amplitudes, low-level talker), (a) 8-digit 
folded binary code, (b) 8-digit binary code, (c) 8-digit Gray code. 

(Hi) No more than one digital error occurs in any given code word. 

(iv) The input signal is speech and is assumed to be Laplace dis- 
tributed. 345 

(/>) The mean value of the input signal corresponds to the mid-range 
of the coder-decoder. 

(vi) Logarithmic companding 12 with /x = 50, 100, and 200. 

(vii) The effects of overload have been ignored. 
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Fig. 9 — Probability of error — /i = 100 log compandor symmetric for nega- 
tive amplitudes average talker, (a) 8-digit folded binary code (speaker level 25 db 
BFLSW). (b) 8-digit binary code (speaker level 30 db BFLSW). (c) 8-digit 
Gray code (speaker level 25 db BFLSW). 
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Fig. 10 — Mean square error vs digit — 2-volt peak-to-peak signal — (a) 
8-digit codes (speaker level 15 db BFLSW — louder talker), (b) 8-digit codes 
(speaker level 30 db BFLSW — average talker), (c) 8-digit codes (speaker level 
50 db BFLSW — low-level talker). 

3.2 Procedure 

3.2.1 Mean Square Error 

A block diagram of the system assumed for the calculations per- 
formed is given in Fig. 1. The results were calculated using a probability 
of error of 10 -6 . This error rate is sufficiently low to justify assumption 
(Hi) above. The probability of a given word occurring on the line depends 
on the compression factor of the coder and the input statistics. The 
mean square error caused by digital line errors for the assumptions 
made is* 



7? 2 = 2-i 2-/ € KJ PkPjLi 



J=\ K M 



1) 



* Formulations have been presented in the literature, different from the 
above, which may be more suitable for determining an optimal nonredundant 
code. 8 - 7 
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where 

t K j is the amplitude of the error for an error in digit J of word K 

P K is the probability of a particular code word occurring 

Pj is the probability of a particular digit of a given word being 

chosen* 
P e is the probability of a digital error (error rate) 
L is the number of digits in the coder 
M = 2 (L ~ l) 
Using the assumptions listed above, we reduce (1) to 

? = P. Z v/Pj (2) 

where rjj 2 is the mean square error for a given digit and is defined in (3) 
as 

^ = 2E^X- (3) 

The calculation of the values given in (2) and (3) are straightforward 
but quite tedious in the case of a compandored system with varying 
speech levels. The problem was programmed for the IBM 7094 computer 
and the results of the computation for (2) and (3) presented in Figs. 4, 
5. and 10(a), (b), and (c). 

3.2.2 Word Probability] 

The probability of a given word K occurring for a given speaker level 
was calculated directly from (4) 2 - 3 ' 4 - 5 



r B K+i 
= C / exp (-2C\x\)dx (4) 



'Br + 1 
'"K 

where 

RMS VALUE FULL LOAD SINE WAVE E 



C = 



RMS VALUE SIGNAL y/2t 



and the B K are normalized coder levels for a specified compression char- 
acteristic. Normalization is with respect to full load. 

The value of P K given above is for a zero mean input signal. When the 
input signal and the coder do not have a zero mean, a shift can be per- 
formed easily. Using the assumptions given, the following result for P K 



* Pj = 1/L for the case considered. 

t An analytic approach is presented in Appendix A. 
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was obtained, 



P K = ffexp (-2CB K ) - exp (-2CB K+ i)}. 



(5) 



The word probability for speaker levels of 20 db, 30 db, 40 db, and 50 
db BFLSW have been plotted in Fig. 3 for n = 50, 100, and 200. 



3.2.3 Line Errors 

The introduction of line errors was accomplished by complementing 
a given digit in a specific code word and determining the new word ob- 
tained. The results obtained obviously depend on the code being con- 
sidered. A transition matrix was produced for each particular code 
considered, i.e., Gray, binary, and folded binary. 

The rows of the transition matrix represent the originally transmitted 
code word and the columns the digit in error from most to least signifi- 
cant. The entry in column J and row K is the new word R produced by 
an error in digit ./ of word K. The transition matrix for a 3-digit binary 
code is given in Table I. This matrix is then used to determine the error 
caused by a specified digital error. The operation required to determine 
e KJ is decoding. If the output voltages corresponding to the code words 
K and R are denoted by D K and D R respectively, we obtain the follow- 
ing expression for e KJ , 

e KJ = D R - D K . (6) 

One should note that the digit in error was specified as digit J in the 
above example. 

Table I — Transition Matrix for 3-Digit Binary Code 











DigitinError 




[_ 


Word -in -Error 


1 o 


1 2 3 


a 
.2 

■4^ 


000 










4 2 1 

5 3 
3 
7 1 2 

5 

1 7 4 

2 4 7 

3 5 G 




— 
- 

a. 


| 001 


| 1 


X 

— 


| 010 


1 2 


~ 


I on 


1 3 


— 
z 


| 100 


| 4 


>> 

— 


| 101 


5 


a 

s 


| 110 

1 111 


1 6 




1 7 















2240 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1965 

The error matrix having elements e K j was then used to calculate the 
mean square error for a given digit, the mean square error for a given 
error rate with and without the "click" component and the distribu- 
tion of error amplitudes for each code considered. Criteria other than 
mean square can be easily adapted to the computation. 

IV. REMARKS AND CONCLUSIONS 

The results obtained indicate the folded binary code yields better 
click performance for a given probability of digital error when compared 
to Gray and binary transmission. The effects of clicks should be virtually 
nonexistent at low talker levels in the folded binary code when compared 
to binary and Gray code. The question arises, however, concerning the 
subjective effect of the line error distortion when compared to quantizing 
error since the clicks at low levels have been suppressed. A 3-db reduc- 
tion in S/N will occur for an 8-digit, /x = 100, 10~ 6 error rate system 
transmitting folded binary for speaker levels between 30 db and 50 db 
BFLSW. The system performance will also be degraded by 3 db if line 
error distortion is subjectively equivalent to quantizing error. 

Further work in evaluating the subjective effect of clicks and noise 
occurring simultaneously is required. 

V. ACKNOWLEDGMENT 

The author wishes to acknowledge the invaluable advice and aid on 
programming given by Miss E. G. Cheathem. 

APPENDIX 

Calculation of Density Function and Word Probability After Compression 

A.l Word Probability 

An analytic expression for the probability density function at the out- 
put of a logarithmic compressor is derived. The density function is de- 
rived for a speech input which can be approximated by a Laplacian 
distribution . 2 - 3 ' 46 The density function used is given in (7) as 

p(ed = C exp (-2C | e t |), (7) 

where, e,- is the input signal and C is denned by 

RMS VALUE OF FULL- LOAD SINE WAVE 



C = 



RMS VALUE OF SIGNAL 
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Normalizing the result to full-load sine wave of unit amplitude results 
in a value \/y/2<i} for C. 

t : ; The characteristic of a logarithmic compressor normalized to a full- 
load unit sine wave is given by 2 

ioga+q) for 0ifia51 ( 8) 

log (1 + m) 
and 

logji^q) for _ u , |0 . (8b) 

log (1 + n) 

The output density function can be determined directly by transform- 
ing the input, i.e., 



dei 

di<o 



(9) 



p(i'o) = p(e,-) 
Performing the operation denned in (9) yields, 

p( „ ) = cio e d+») (1 + m)( . exp r_?c l(1 + m) ,, _ lt i 

M L m J (10) 

for n > and ^o ^ 

andp(-fo) = p(y ). The expression given in (10) is the probability density 
function of the compressor output. 

Determining the word probability for a given code word A' involves 
integrating the derived density between the lower and upper transition 
values for which this code word is emitted. For the case of a large num- 
ber of digits (fine quantizing), an approximate expression can be de- 
rived. The step size for an n digit quantizer having a peak-to-peak 
amplitude of 2 is l/2 ( "-°. The value of the word probability P(K) for 
word K is given by; 

p(g) = cjogji + M) r-- (1 + M) .. 

H • / (K-l)/2 n_1 



, (ID 

dVa 



.exp[-—{(l +/*r- 1} 

for K ^ 1 and v ^ 0. 
We can approximate this result for the case of fine step sizes by 
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P(K) m Clog(l+») (1 + j e «-„«»-. 

2" V 

for K ^ 1 and v ^ 0, 

or 

P(K) S Cl0g o (1 1 + /t) (1 + M ) ^ (2 "" I)/2n 
2" > 

•exp[-^[(l + M ) e (2K -'" 2n -l|]. 

For the case of n = 8, p = 100, # = 1, C = 100 (40 db BFLSW) 

P(l) ^ ^ (1.002)<f 2(0 002) ^ ^ (1.002) S 0.036. (14) 
12o 128 

This agrees favorably with the value computed from the exact equation 
for this case of 0.0351 . 

A. 2 Determination of the Most Probable Code 

The most probable code level will correspond approximately to the 
peak(s) of the density function when the quantizing steps are small. For 
this condition we obtain, 

g*> = . Cln'd+M) (1 + „)... exp [_ 20 

dv n L M 

•[(1 + M)"° m - 1]] X |l - y (1 + M)^} (15) 

for v^ ^ 
where v 0m is the peak value of v . Hence we have, 

A - u + ,)- 



or 



log (1 + n) 
The results obtained using (16) compare favorably with those deter- 
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Table II — Comparison of Computed and Analytically 
Obtained Values 
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mined on the computer as indicated in Table II. The value given for 
v 0m is the decimal fraction of full load code levels. 
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